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Abstract.  The  main  aim  of  this  article  is  to  investigate  some  new  families  of  g-special 
polynomials  and  to  study  their  properties  using  different  approaches.  The  2-iterated 
q-Appell  polynomial  family  is  introduced  and  some  properties  of  these  polynomials  are 
considered  under  q-umbral  calculus  methods.  Some  2-iterated  q-Appell  and  hybrid  q- 
special  polynomials  are  studied  as  members  of  this  family.  The  numbers  related  to  these 
polynomials  are  obtained.  The  graphical  representation  of  the  2-iterated  and  hybrid  q- 
Appell  polynomials  is  presented.  The  zeros  of  these  polynomials  are  investigated  for  some 
values  of  index  n  using  numerical  computation.  The  approximate  solutions  of  the  real 
zeros  of  these  polynomials  are  also  considered.  The  determinant  forms  for  the  2-iterated 
(j-Appell  family  and  for  the  2-iterated  and  hybrid  g-special  polynomials  are  established 
using  linear  algebraic  approach. 
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1.  Introduction  and  preliminaries 

The  area  of  g-calculus  has  in  the  last  twenty  years  served  as  a  bridge  between 
mathematics  and  physics.  Recently,  there  has  been  a  significant  increase  of 
activities  in  the  area  of  g-calculus  due  to  its  applications  in  various  fields  such 
as  mechanics,  mathematics  and  physics. 

The  definitions  and  notations  of  g-calculus  reviewed  here  are  taken  from  [3]. 
The  g- analogue  of  the  shifted  factorial  {a)n  is  defined  by 
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n— 1 


(a;  q)o  =  1,  (a;  q)n  ^  ^ 


m=0 


The  g-analogues  of  a  complex  number  a  and  of  the  factorial  function  are 
defined  by 

1-0“ 

[a]g  =  ^gC-{1};  a  G  C, 


w,!  =  n  H,  =  iii,[2i,  ■  ■  ■  w,  = 

m=l  ' 

q  ^  1;  n  G  N,  [0]^!  =  1,  (7  G  C;  0  <  (7  <  1. 
The  Gauss  g-binomial  coefficient  [^]  ^  is  defined  by 


nl  ^  [n]ql  ^  {q;q)n 
k\q  [k]ql[n-k]ql  {q]q)k{q;q)n-k 

The  (7-exponential  function  is  defined  as: 


A:  =  0, 1, . . . ,  n. 


00 


n=0 


0  <  |g|  <  1. 


The  (7-derivative  Dqf  of  a  function  /  at  a  point  0  7^  z  G  C  is  defined  as: 

qz  —  z 

The  (7-analogue  of  Taylor  series  expansion  of  an  arbitrary  function  f(z)  for 
0  <  g  <  1  is  defined  as: 


00 


/(^)  =  E 

n=0 


(1-gr 

{q;q)n 


D^,f{a){z 


a) 


n 


where  D^f{a)  is  the  (7-derivative  of  the  function  /  at  point  a. 

In  1985,  Roman  proposed  an  approach  similar  to  the  umbral  approach  under 
the  area  of  nonclassical  umbral  calculus  which  is  called  (7- umbral  calculus  [11,  9]. 
Let  P  be  the  algebra  of  polynomials  in  the  single  variable  x  over  the  field  of 
complex  numbers.  Let  P*  be  the  vector  space  of  all  linear  functionals  on  P.  Let 
{L  I  p{x))  be  the  action  of  a  linear  functional  L  on  a  polynomial  p{x).  Let  5 
denote  the  algebra  of  formal  power  series 


00 

m  =  E 


k=0 
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This  kind  of  algebra  is  called  an  umbral  algebra.  Each  f{t)  £  ^  defines  a 
linear  functional  on  P  given  by  {f{t)\x^)  =  for  all  A:  >  0.  In  particular, 
=  n\5n,k  for  all  n,k  >  0,  where  Sn,k  is  the  Kronecker  delta.  Let  Cn  be 
a  fixed  sequence  of  nonzero  constants.  Then  the  algebra  5  of  all  formal  power 
series  in  t  is  isomorphic  onto  the  vector  space  P*  by  setting 


If 


(A  I  a;”)  =  Cn5n,k 

OO 

fit)  = 


OO  ^k 


k=0 


Ck 


then 


{fit)\x^)  =  Ofc. 

The  g-umbral  calculus  is  defined  by  setting 


C-n.  — 


(1  -  qy 


Then 


Cn 

Cn— 1 


The  functional  eq{yt)  satisfies 


1-q^ 

1-q' 


ie,iyt)\xl  =  =  y" 


Then 


{eqiyt)\pix))  =p{y), 

for  all  polynomials  p{x)  G  P.  Also,  the  operator  eq{yt)  satisfies 


eq{yt)x^  = 


h 

k=o  1 


n 


/  x^-y 


Let  f(t),  g{t)  G  5^.  Then  we  have 


ifit)git)\pix))  =  {fit)  I  git)piy)). 


In  1967,  Al-Salaam  [1]  introduced  the  family  of  g-Appell  polynomials 
{An,q(x)}n>o  and  studied  some  of  its  properties.  The  n-degree  polynomials 
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-^n,q{x)  are  called  (7-Appell  if  and  only  if  these  are  defined  by  means  of  the  fol¬ 
lowing  generating  function  [1]: 


Aq(t)eq(xt)  =  V  A„,g(x)-pp-,  0  <  g  <  1,  (1) 

where 

^q{t)  ■=  ^n,q  r  i  |  ;  ^0,q  =  Ij  ^q{t)  7^  0.  (2) 

n=0 

It  is  to  be  noted  that  Aq{t)  is  an  analytic  function  at  t  =  0  and  An^q  :=  An^q{0). 
Also,  there  exists  a  sequence  of  numbers  {An^q}n>o,  such  that  the  polynomials 
An,q{x)  satisfy  the  following  relation  [1]; 


An,q{x)  —  An,q  + 


n 

J 


An-l^q  X  + 


n 


An-2, q  X^  +  ■  ■  ■  +  Aq  q  x”,  n  —  0,  1,  2,  .  .  .  . 


The  g-Appell  polynomials  An,q{x)  are  characterized  by  Roman  [10]  under 
g-umbral  calculus.  Let  An,q{x)  be  g-Appell  polynomials  for  g{t).  Then 


n=0 


^n,q{y) 


N<?! 


To 


9q{tW‘ 


The  polynomials  An,q{x)  are  the  (7-Appell  polynomials  for  gq{t)  if  and  only  if 


i. 

1  ^  -^n 

ii. 

An,q{x)  =  {gq{t))~^X^,  (3) 

in. 


eq{yt)An,q{x)  = 


fc=0  L  -'9 


An—k,q{x)y 


IV. 
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Al-Salaam  showed  that  the  class  of  all  g-Appell  polynomials  is  a  maximal  com¬ 
mutative  subgroup  of  the  group  of  all  polynomial  sets,  i.e.  the  class  of  all  g-Appell 
sequences  is  closed  under  the  operation  of  g-umbral  composition  of  polynomial 

n  n 

sequences.  If  An^q{x)  =  0'n,k-,q  and  Bn,q{x)  =  ^  bn,k;q  x^  are  sequences  of 

k=0  k=0 

g-polynomials,  then  the  g-umbral  composition  of  An^q{x)  with  Bn,q{x)  is  defined 
to  be  the  sequence 

n 

{An,q  O  Bq){pc)  =  ^  ^  Q'n,k\q  B}^^q{x)  =  ^  ^  0‘n,k\qbk,l\qX  ■ 

k=0  0<k<l<n 

Under  this  operation,  the  set  of  all  g-Appell  sequences  is  an  abelian  group 
and  it  can  be  seen  by  considering  the  fact  that  every  g-Appell  sequence  is  of  the 
form 

and  that  the  umbral  composition  of  g-Appell  sequences  corresponds  to  multipli¬ 
cation  of  these  formal  g-power  series  in  the  operator  Dq. 

Based  on  appropriate  selection  for  the  function  Ag{t)  or  different 

members  belonging  to  the  family  of  g-Appell  polynomials  can  be  obtained. 
These  members  are  listed  in  Table  1. 


Table  1.  Some  members  belonging  to  the  g-Appell  family 
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The  article  is  organized  as  follows.  In  Section  2,  the  2-iterated  g-Appell 
polynomials  are  introduced  by  means  of  generating  function  and  series  definition. 
Some  properties  of  these  polynomials  are  considered  under  g-umbral  calculus. 
The  members  belonging  to  the  2-iterated  (7-Appell  family  and  some  hybrid  q- 
special  polynomials  and  related  numbers  are  also  explored.  In  Section  3,  the 
zeros  of  the  2-iterated  g-Appell  and  hybrid  g-Appell  polynomials  are  investigated 
using  numerical  computation  and  their  graphs  are  drawn  for  suitable  values  of 
indices.  In  Section  4,  the  determinant  forms  of  the  2-iterated  and  hybrid  (7-Appell 
polynomials  are  established. 


2.  2-iterated  g-Appell  polynomials 

In  order  to  introduce  the  2IgAP,  two  different  sets  of  g-Appell  polynomials 
g(x)  and  are  considered.  Thus,  from  dehnitions  (1)  and  (2),  it  follows 

that 

A\{t)eq{xt)  =  0<q<l, 

n=o 


where 


(4) 


j-n 

A^it)  :=  ^  :=  A^g(O);  A^^  =  1;  A^(t)  /  0 

n=o 


(5) 


and 


where 


A\\t)eq{xt)  =  0  <  g  <  1, 


n=0 


(6) 


^  fTl 

Aq^it)  :=  ^n!g(0);  =  1;  /  0- 

n=0 


2.1.  Generating  function  and  other  properties 

The  generating  function  for  the  2IgAP  is  derived  by  using  a  different  approach 
based  on  replacement  techniques.  For  this,  the  following  result  is  proved; 
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\2] 

Theorem  1.  The  polynomials  Anjq{x)  are  the  2-iterated  q-Appell  polynomials  for 
Aq{t)  and  A^^  ft)  (or  for  (t)  and  g^^ft))  if  and  only  if 


Gq{x,t)  :=  Al{t)Al\t)eq{xt)  =  A]^}^{x) ,  0  <  q  <  I,  (7) 

n=o 

or,  equivalently, 

1 

g({t)g(i ^  °  ^  ^  ^ 


Proof.  Expanding  the  (/-exponential  function  eq{xt)  in  the  l.h.s.  of  equation 
(4)  and  then  replacing  the  powers  of  x,  i.e.  x^ ,  x^,  . . . ,  by  the  correspond¬ 

ing  polynomials  ^^^(x),  A\^^q{x),  . . . ,  Al^q{x)  in  the  l.h.s.  and  replacing  x  by  the 
polynomial  AY^{x)  in  the  r.h.s.  of  the  resultant  equation,  it  follows  that 


r  t  in 

1  +  ^lfg(^)'jYr7  +  +  •  •  •  +  Mi,q{^)YY~\  +  '  '  ' 

L  [^q-  [n\q. 


n 


I 


n=0 

Summing  up  the  series  in  l.h.s.  and  then  using  equation  (6)  and  denoting  the 

\2] 

resultant  2IgAP  in  the  r.h.s.  by  An)q{x),  that  is 

^n}qix)  =  A(^g{A{^^g{x)},  (9) 

assertion  (7)  is  proved.  In  view  of  equations  (1)  and  (3),  assertion  (8)  follows.  ◄ 

Next,  the  series  definition  for  the  2-iterated  (/-Appell  polynomials  is  obtained 
by  proving  the  following  result: 

\2] 

Theorem  2.  The  polynomials  Anjq{x)  are  the  2-iterated  q-Appell  polynomials  for 
A^  q(x)  and  A(^g(x)  if  and  only  if 

^n,U*)  = 


n  r  -I 

n 


k=0 


k 


An-k,q{x) 


(10) 


Proof  Using  equations  (5)  and  (6)  in  the  l.h.s.  of  generating  function  (7)  and 
then  using  Cauchy-product  rule  in  the  l.h.s.  of  the  resultant  equation,  we  obtain 


EE 

n=0  k=0 


n 


K,q  K-k,q{x) 


I  I  I  A^-n,q 

^  n=0 


X] 


N<?! 


(11) 
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Equating  the  coefficients  of  same  powers  of  t  in  both  sides  of  equation  (11), 
assertion  (10)  follows.  ◄ 

Further,  as  consequences  of  some  results  mentioned  in  Section  1,  the  following 
theorems  are  deduced; 

[21 

Theorem  3.  The  polynomials  An'q{x)  are  the  2-iterated  q-Appell  polynomials  if 
and  only  if 


{eq{yt)  I  Ag'g(x))  =  I  =  ^n]q{y)- 

\2] 

Theorem  4.  The  polynomials  Anlq{x)  are  the  2-iterated  q-Appell  polynomials  for 
gUt)  and  g^^it),  if  and  only  if 


liyi)  = 


(eqiyt)  I  A^^^{x))  J  JI  f,  _  \,q{y)  J  (.-sj!  f.-s.k 


.[2] 


k=0 


[kW- 


-9q{t)9q\t)t^  =  -T§n-9q{t)9q\t)t 

fc=o 


f2l 

Theorem  5.  The  polynomials  An;q{x)  are  the  2-iterated  q-Appell  polynomials  for 
gq{t)  and  gq\f),  if  and  only  if 

2i£Ua;)  =  {9!,{t)gl^{t))-^  x^. 


\2] 

Theorem  6.  The  polynomials  An'q{x)  are  the  2-iterated  q-Appell  polynomials,  if 
and  only  if 


t  ^n]q{x)  =  [n]q  7ljfli_q(x). 


[21 

Theorem  7.  The  polynomials  An;q{x)  are  the  2-iterated  q-Appell  polynomials,  if 
and  only  if 


eq{yt)A^^\{x)  =  ^ 


fc=o 


f2l 

Theorem  8.  The  polynomials  Anjq{x)  are  the  2-iterated  q-Appell  polynomials  for 
A(^  g(x)  and  Al^q{x)  if  and  only  if 
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2.2.  Particular  members 

By  making  suitable  selections  for  the  functions  and  Ag^{t)  in  equations 

(7)  and  (10),  the  generating  function  and  series  definition  for  the  corresponding 
member  belonging  to  the  2-iterated  g-Appell  family  can  be  obtained.  These 
resultant  2-iterated  (7-Appell  polynomials  along  with  their  notations,  names, 
generating  functions  and  series  definitions  are  given  in  Table  2. 

Table  2.  Some  members  belonging  to  the  2-Iterated  g-Appell  family 


S.  ^g(i)  =  Resultant  Generating  function  Series  definition 

No. _ 2IgAP _ 


I. 


g(t)-l 


=  BS}g(x)  = 

2-iterated  =  (a:)  [fc]  ^  S„  _ (a;) 

q- Bernoulli 
polynomials 

(2IgBP) _ 


II. 


III. 


( 


g(t)  +  l 


) 


2-iterated  =  (a:)  EIE  =  ijl  Bfc,,  fc ,,  (a:) 

ra=0  '-"■'1-  k=0 

g-Euler 

polynomials 

(2IgEP) _ 


E  (-1)” 

n  =  0 


rl(>l-l)  _ 


r[21 


(x)-- 


2-iterated 

g-Hermite 

polynomials 

(2IgHP) 


E  (-1)"?' 


Xeq(xt)  =  y;  Hl^}g(x) 
n  =  0 


BS}g(x) 

=  E  [jl  Bf.  (x) 


The  combinations  of  any  two  different  members  of  the  g-Appell  family  in 
the  2-iterated  (7-Appell  family  yields  a  new  hybrid  (7-special  polynomial.  Thus, 
by  making  suitable  selections  for  the  functions  Ag(t)  and  Ag^{t)  in  equations 
(7)  and  (10),  the  generating  function  and  series  definition  of  the  resultant 
hybrid  (/-special  polynomials  can  be  obtained.  The  possible  combinations  of  the 
g-Bernoulli,  g-Euler  and  g-Hermite  polynomials  (Table  1  (I-III))  are  considered. 
The  resultant  hybrid  (/-special  polynomials  are  given  in  Table  3. 


Remark  1.  Note  that 


—  EBn^qix)  t  BHn,qix)  —  H  Bn^q{x^ ,  —  JJ  E^iqi^X^ , 

where  EBn,q{x),  HBn,q{x)  and  HEn,q{x)  are  the  q- Euler- Bernoulli  polynomials 
(qEBP),  q-Hermite- Bernoulli  polynomials  (qHBP)  and  q-Hermite- Euler  polyno¬ 
mials  (qHEP).  This  confirms  the  faet  that  the  set  of  all  q-Appell  sequenees  is 
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closed  under  the  operation  of  q-umbral  eomposition  of  polynomial  sequences  and 
forms  an  abelian  group. 

Table  3.  Some  hybrid  g-special  polynomials 


s. 

No. 

I. 


Hybrid  Generating  function  Series  definition 

polynomial 

BEn,q(x)  ~  (e,(t)-l)(eg(t)  +  l)  = 

q-Bernoulli-  =  2  BEn,q(x)j:^  =  E  [  j]  ^  Bfc,,  Bn  -  fc,,  (a;) 

Euler 

polynomials 

(gBEP) _ 


II. 


(  eq(t)-l  )  ’ 


BHn,q(x)~ 

g-Bernoulli- 

Hermite 

polynomials 

(gBHP) 


(^qd)-O 

i2n 


[2n],!! 

=  E  BH„,g(x) 
n=0 


j(xt) 


BHn,q(x)  - 

=  E  [fcl  g(x) 

k  =  0  '■  ■'  1 


III. 


(i^);  eG„.,W:  = 

oo 

E  (-1)"  -  O-Euler- 

n  =  o'  [2«],!! 

Hermite 


polynomials 

(gEHP) 


[2n],_ 

=  £  Ef^n,q(x) 

n=0 


-eq(a:£) 


£;Hn,q(a:)  = 


2.3.  Related  numbers 

f2l 

The  numbers  related  to  the  members  of  the  2IqAP  An^q{x)  given  in  Tables 

f2l 

2  and  3  are  obtained.  By  taking  re  =  0  in  series  definitions  of  2IqBP 
2IqEP  En}q{x)  and  2I(7HP  Hn}q{x)  given  in  Table  2  (I-III)  and  using  notations 
from  Table  1,  the  2-iterated  q-Bernoulli,  q-Euler  and  q-Hermite  numbers  are 
obtained.  These  numbers  are  given  in  Table  4. 

Table  4.  2-Iterated  q-numbers 


S.No.  Notation  and  name  of  the  2-iterated  g-number  Series  definition 


I. 

Bf]q  :=  bP>,(0) 

2-iterated  g-Bernoulli  numbers  (2IgBN) 

si?,'. 

II 

\l]^Bk,qBn-k,q 

II. 

Ef:]q  -  bPi,(o) 

2-iterated  g-Euler  numbers  (2IgEN) 

II 

III. 

Hf}q  ■.=  B|f',(0) 

2-iterated  g-Hermite  numbers  (2IgHN) 

II 

[fc]  ^^k,q^n-k,q 
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Similarly,  the  numbers  corresponding  to  the  hybrid  g-special  polynomials 
given  in  Table  3  can  be  obtained.  Taking  x  =  0  in  series  definitions  of  the 
gBEP  BEn,q{x),  gBHP  BHn,q{x)  and  (7EHP  EHn,q{x)  (Table  3  (I-III))  and 
using  notations  from  Table  1,  the  (7-Bernoulli-Euler,  g-Bernoulli-Hermite  and 
(7-Euler-Hermite  numbers  are  obtained.  These  numbers  are  given  in  Table  5. 

Table  5.  Hybrid  g-numbers 


S.No.  Notation  and  name  of  the  mixed  type  g-numbers  Series  definition 


I. 

BEn,q  ~  B-®n,q(0) 

g-Bernoulli-Euler  numbers  (gBEN) 

BEn,q  ^  S 

fc=0 

[J] 

Ek,qBn-k,q 

II. 

g-Bernoulli-Hermite  numbers  (gBHN) 

n 

Bfin,q  = 

fc=0 

[2] 

Ek.q  Efx  —  k.q 

'  1 

III. 

EHn,q  ~  EJfn,g(0) 

g-Euler-Hermite  numbers  (gEHN) 

II 

[fe] 

^  Ek,qEfi  —  k,q 

Remark  2.  Also  note  that 

sEn^q  —  EBn,q  EBn,q{0)',  sHn^q  —  uBn^q  •  = 

:=  HBn,q{0)-,  EHn,q  =  HEn,q  ■=  HEn,q{0),  (12) 

where  EBn,q,  nBn^q  and  HEn,q  are  the  q- Euler- Bernoulli  numbers  (qEBN),  q- 
Hermite- Bernoulli  numbers  (qHBN)  and  q-Hermite- Euler  numbers  (qHEN),  re¬ 
spectively. 

In  the  next  section,  the  shapes  of  the  2I(7BP  Bn}q{x),  2IgEP  En}q{x),  2IgHP 
Hn}q{x),  (7BEP  BEn,q{x),  (7BHP  BHn,q{x)  and  (/EBP  EHn,q{x)  are  displayed.  The 
zeros  of  these  polynomials  are  also  investigated  using  numerical  computation. 

3.  Zeros  and  approximate  solutions 

There  has  been  increasing  interest  in  solving  mathematical  problems  with  the 
aid  of  computers.  The  numerical  investigation  of  the  zeros  of  some  (/-polynomials 
are  considered  in  [12,  13,  14,  7].  First,  we  give  the  shapes  of  the  g-special  poly¬ 
nomials. 

3.1.  Shapes  of  the  2-iterated  and  hybrid  g-Appell  polynomials 

The  plots  of  the  2IgBP  B^^q{x),  2IgEP  E^^q{x),  2IgHP  H^^q{x),  gBEP 
BEn,q{x),  gBHP  BHn,q{x)  and  (/EHP  EHn,q{x)  are  drawn  for  n  =  1,2,  3, 4  and 
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(7  =  ^  (0  <  g  <  1).  This  shows  the  four  plots  combined  into  one  for  each  of  these 
polynomials.  For  this,  the  values  of  the  first  five  Bn^q,  En,q  and  Hn,q  are  re¬ 
quired.  The  values  of  first  five  B^^q,  E^^q  and  Hn,q  [4,  6]  are  mentioned  in  Table  6. 


Table  6.  Values  of  first  five  Bn^q,  E^^q  and  Hn,q 


n 

0 

1 

2 

3 

4 

Sn,, 

1 

-(1  +  9)”" 

P([3],!)”^ 

(1 

-9)P([2],)-h[4],)-l 

9hl  -  9"  -  292  -  9W  9®)([2]^[3],[51,)-1 

1 

1 

■j(  — 1  +  q) 

-1  +  29  +  292  -q3) 

4(9-  1)[3]9!(9^  -49  +  1) 

1 

0 

-1 

0 

[3)99" 

The  expressions  of  the  first  five  Bn,q{x),  En,q{x)  and  Hn,q{x)  are  obtained 
by  making  use  of  the  values  of  the  first  five  Bn,q,  En^q  and  Hn^q  in  the  series 
definitions  given  in  Table  1  (I-III).  The  expressions  of  first  five  Bn,q{x),  En,q{x) 
and  Hn,q{x)  are  mentioned  in  Table  7. 


Table  7.  Expressions  of  first  five  Bn,q{x),  En,q{x)  and  Hn,q{x) 


n 

0 

1 

2 

3 

4 

Bn,q(x) 

1 

X  — 

1 

1  +  q 

2  _  PI9  , 

l  +  q  4- 

+  [3lq[2]. 

3  [3]q^^  1  q'^x  1  (l-g)g^ 

1  +  9  +  T2I7  121 9 1+9 

4  [4]q  3  ,  [4]qg^  2  , 

1=  i+l^  +  + 

+  9^(1  -  9^  -  29^- 
-9W9®)([2]^[3]q[5]q)-l 

En,q(x) 

1 

X  — 

1 

2 

^2  _  + 

+  1  (  — 1  +  9) 

+  ^(-1  +  29  +  292  -q3) 

^4  1+99,3^  [4191319(9-1)^2^ 

14]9(-1  +  29+292-93) 

g  0;  + 

,  (q-l)[3]q!(q2_4q+l) 

Ifi 

Hn,q(x) 

1 

X 

9,2  _  1 

-  [3],  X 

x‘^  -  (1  +  9^)  [3],  x'^  +  [3]^  9 

By  making  appropriate  substitutions  from  Tables  6  and  7  in  the  series 
definitions  given  in  Tables  2  and  3,  the  expressions  of  the  first  five  2I(7BP 
Bl^}q{x),  2IgEP  El^}q{x),  2IgHP  H^q{x),  gBEP  BEn,q{x),  g'BHP  BHn,q{x)  and 
gEHP  EHn,q{x)  for  (7  =  ^  are  obtained.  These  expressions  are  given  in  Table  8. 

With  the  help  of  Matlab  and  by  using  the  expressions  for  the  first  five 

^ip/2(®)’  BEn^i/2ix),  bH^, i/2ix)  and  eH^, i/2ix)  from 
Table  8  for  n  =  1,  2, 3, 4,  the  following  graphs  are  drawn: 


32 


M.  Riyasat,  S.  Khan,  N.I.  Mahmudov 


Shape  of 


Table  8.  Expressions  of  first  five 

EHn^l/2{x) 


Shape  of 


Figure  3.2 


(®)’  ^n,l/2(®)’  ^S/2(®)’  B^n,l/2(a;), 


n 

0 

1 

2 

3 

4 

1 

X  — 

4 

3 

-  2a:  +  f 

a:3  - 

Ia:2  ,  3^  _  _8_ 

3X  -r  2-^  45 

- 

5 „3  ,  45 „2 

2  ^24 

8  „  221 

24  7812 

1 

X  — 

1 

-  ^ 

2  16 

a:^  - 

Ia;2  ,  ,  _5_ 

4  ~  32  ~  16 

a:^  - 

15  ^3 
8  ^ 

^  128 

+  125  ,  71 

^  1024  ^  1024 

<\/2(-) 

1 

X 

-  2 

- 

^x 

2  ^ 

a:^  - 

35  ^2 
8  ^ 

1  56 
■'■  16 

,,1/2(2:) 

1 

X  — 

7 

fi 

_  3^,  , 

a:3  - 

49  „2  ,  79^  ,  379 

24  X  -r  X  -r  28sn 

a:*  - 

M^3 

16 

M5  ^2 
384 

461  ,,  402305 

1536  9999360 

.,1/2(2:) 

1 

X  — 

3 

? 

-  X  - 

21 

a:3  - 

7„2  19  „  ,  107 

6  1 2  ~  90 

- 

^x^  - 
4^ 

100^2  ,  107^  ,  14554 

48  ^  48  ^  62496 

.,1/2(2:) 

1 

X  — 

1 

2 

-  ^X  -  ^ 

4-^  8 

a:3  - 

7„2  63  „  ,  59 

8  32  ^  64 

a:‘‘  - 

15  ^3 

16 

315  ^2 
128 

4_  885  „  4_  791 
^  512-^  ^  1024 

Shape  of  Hra^^(x)  Shape  of 


Figure  3.3 


Figure  3.4 
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Figure  3.6 


3.2.  Zeros  and  approximate  solutions  of  the  2-iterated  and  hybrid 
g-Appell  polynomials 

We  find  the  real  and  complex  zeros  of  the  polynomials  given  in  Table  8.  The 
manual  computation  of  these  zeros  is  too  complicated,  therefore,  we  use  Matlab 
to  investigate  these  zeros.  The  investigation  in  this  direction  will  lead  to  a  new 
approach  employing  numerical  methods  in  the  field  of  these  g-special  polynomials. 

The  real  zeros  of  BEn, i/2{x),  i/2{x)  and 

EHn^i/2ix)  are  computed  by  using  Matlab.  These  zeros  are  given  in  Table  9. 

Again,  with  the  help  of  Matlab,  we  find  the  complex  zeros  of  these  polyno¬ 
mials.  These  complex  zeros  are  given  in  Table  10. 


Table  9.  Real  zeros  of  i?F„^i/2(a;),  i?R„^i/2(x), 

1/2(2:) 


Degree  n 

<’l/2C) 

^n,'l/2C) 

BEn,l/2(^) 

1/2(2:) 

1 

1.3333 

1.0000 

0.0000 

1.1667 

0.6667 

0.5000 

2 

0.6220,  1.3780 

-0.0406,  1.5406 

1.4142,  -1.4142 

0.2411,  1.2589 

1.5746,  -0.5746 

1.5000,  -0.7500 

3 

0.1522,  0.9446, 
1.2365 

-3.2768,  -0.2642, 
2.3743 

0.0000,  1.8708 
-1.8708 

-0.1213,  0.7910, 
1.3719 

-1.1392,  1.6874 
0.6184 

-1.2626,  1.7108, 
0.4268 

4 

-0.0617,  0.3823 

0.5479,  1.6488 

-1.8224,  -1.0266 
1.8224,  1.0266 

-0.2476,  2.6664 

-1.3694,  1.6381, 
1.0777,  -0.0963 

-1.4071,  1.6153, 
1.0523,  -0.3230 

Remark  3.  From  Tables  9  and  10,  the  following  general  relation  is  observed. 
The  number  of  real  zeros  lying  on  the  real  plane  Im  (x)  =  0,  i.e., 
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Real  zeros  of  A]^^g{x)  =  n  —  Complex  zeros  of  A]^^g{x), 
where  n  is  the  degree  of  the  polynomial. 

Table  10.  Complex  zeros  of  -^id/2(^)’  s-®n, 1/2(2^)? 

bHu, i/2{x)-,  eHti, i/2{x) 


Degree  n 

BBn,l/2(^) 

B-ffn, 1/2(2:) 

1 

- 

- 

- 

- 

- 

- 

2 

— 

— 

- 

— 

- 

— 

3 

- 

- 

- 

- 

- 

- 

4 

1.0897  -1-  0.1112i, 
1.0897  -  0.1112i 

-0.1609  +  0.225f, 
-0.1609  -  0.225i 

~ 

-0.1157  -1-  0.2181i, 
-0.1157  -  0.2181i 

In  order  to  make  the  above  discussion  more  clear,  we  draw  the  combined 
graphs  of  shapes  and  zeros  of  the  polynomials  i/2(®)’  '^i^i/2(^)’ 

BEn,i/2{x),  bH^, i/2{x)  and  eHu, i/2ix)  for  n  =  4. 

180 
160 
140 
120 
100 
80 
60 
40 
20 
0 

-20 


Shape  and  zeros  of  of  Shape  and  zeros  of  of 


Figure  3.7  Figure  3.8 


Shape  and  zeros  of  of  Shape  and  zeros  of  of  BE4 


Figure  3.9 


Figure  3.10 
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Shape  and  zeros  of  of 


Re(x) 

Figure  3.11 


Shape  and  zeros  of  of 


Re(x) 

Figure  3.12 


It  should  be  noted  that  in  Figures  3.7,  3.8  and  3.10  out  of  total  two  complex 
zeros  only  one,  with  positive  imaginary  part  is  visible,  due  to  the  absence  of 
negative  imaginary  axis  in  these  graphs. 

f2l 

The  numerical  results  for  approximate  solutions  of  real  zeros  of 

sEn, i/2{x),  sI7n,i/2(a^)  and  i/2{x)  for  (n  =  1,2, 3, 4) 
are  displayed  in  Table  9. 

Also,  we  note  that  the  real  zeros  of  these  polynomials  as  shown  in  Table  9  give 

f2l  [2l 

the  numerical  results  for  the  approximate  solutions  of  =  0,  El^^^i^{x)  = 

f2l 

0’  ^n  i/2(^)  =  S-E^n, 1/2(2^)  =  0,  siEn, 1/2(3^)  =  0  and  EiEn,i/2(x)  =  0  for  n  = 

1,2,  3,4. 

In  the  next  section,  the  determinant  forms  for  the  2-iterated  g'-Appell  poly¬ 
nomials  and  some  hybrid  g-special  polynomials  are  established. 


4.  Determinant  forms 

The  g-Appell  polynomials  are  studied  using  determinant  and  umbral  ap¬ 
proaches,  see  for  example  [5,  6].  Obtaining  determinant  forms  for  the  g'-Appell 
polynomials  and  their  members  is  an  important  aspect  of  such  study.  The 
determinant  forms  can  be  helpful  for  computation  purposes  and  can  also  be 
used  in  finding  the  solutions  of  general  linear  interpolation  problems.  This 

fact  provides  motivation  to  establish  the  determinant  forms  of  the  g-special 

f2l 

polynomials  introduced  in  previous  section.  In  order  to  define  the  2Ig'AP  An;q{x) 
by  means  of  determinant,  the  following  result  is  proved: 
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4.1.  Determinant  form  for  the  2-iterated  g-Appell  polynomials 
Theorem  9.  The  following  determinant  form  for  the  2-iterated  q-Appell  polyno- 


2 

mials  An'q 

(x)  of  degree  n 

holds  true: 

Ao^qix)  = 

1 

■  h0,g  ’ 

1 

Ai,q{x) 

A2,q{x)  ■■ 

AiLl,qix) 

An,q{x) 

Po,q 

/3l,q 

1^2, q 

I3n—l,q 

l^n,q 

An}q{x)  = 

.  (-1)" 

0 

f3o,q 

[  1  ]g  l^n-2,q 

[l]g  l^n-l,q 

■  (A, 9)"+^ 

0 

0 

f3o,q 

[  2  ]q  (^n-3,q 

[2]^  Pn-2,q 

0 

0 

0 

Po,q 

L^],  /5i,. 

(13) 

where  n  =  1, 2, . . .  and  {n  =  0, 1, 2, . . .)  are  the  q-Appell  polynomials  of 

degree  n;  /3o,q  ^  0  and 


/3o,g- 

Pri,q  — 


0,9 


n  =  1,2,...  . 


f2l 

Proof.  Let  An;q{x)  be  a  sequence  of  the  2IgAP  defined  by  equation  (7)  and 
^n,q:  l^n,q,  be  two  nuHiericai  sequences  (the  coefficients  of  g-Tayior’s  series  ex¬ 
pansions  of  functions)  such  that 


"4q(t)  —  Aq  g  -I-  -I-  g  H  h  g  H  ,  n  —  0, 1, ... ;  Aq  ^  /  0, 

(14) 

7^  0, 

(15) 


w 


[!]«' 


t 


[2].’ 


f' 


^qi^)  —  Po,q  +  7:j-rTA,q  +  H - +  -f-j— +  '  '  '  ,  n  —  0,  1,  ...  ;  /3o,q  7^  0, 

l^lq- 


n  L! 


satisfying 

Then,  according  to  the  Cauchy-product  rule,  it  follows  that 


Aq{t)A[{t)  =  1. 


n=0  fc=0  ^  ^ 


Ak,q  Pn—k,q 


W' 
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which  gives 

n 

E 

A:=0 

Consequently,  the  following  holds: 


l^n—k,q  — 


if  n  =  0, 
if  n  >  0. 


n  =  1,2, . . . 


Multiplication  of  equation  (7)  by  Ag{t)  gives 


(16) 


{t)eq{xt)  = 


n=0 


n! 


which  in  view  of  equations  (14),  (15)  and  (6)  gives 


n=0 


oo 


n=0 


oo 


y~!  I^k,q 

k=0 


Again,  multiplication  of  the  series  on  the  r.h.s.  of  equation  (16)  according  to 

Cauchy-product  rule  leads  to  the  following  system  of  infinite  equations  in  the 
f2l 

unknowns  An'q{x)  (n  =  0, 1, . . .): 


< 


^ol(a^)/5o,g  =  1, 

^o!qi^)l3l,q  +  Afj^{x)l3o,q  =  A{^,j{x)  , 

^o!qix)/32,q  +  [3g^qg(a;)/3i,g  +  A^^|^(x)/3o,q  =  A^^^qi^), 

^o!qi^)l^n-l,q  +  ["'i  (ic)/3n-2,g  H  +  ^n,g(®)/5o,(j  =  ^n-l,g(2^)) 

^o!qi^)l3n,q  +  [1]  ^Afj^{x) j3n-l,q  +  •  •  •  +  AL%(x)/3o,q  =  Al/g{x), 


First  equation  of  system  (17),  proves  the  first  part  of  assertion  (13).  Also,  the 
special  form  of  system  ((17)  (lower  triangular)  allows  to  work  out  the  unknowns 
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\2] 

A.njg{x).  Operating  with  the  first  n+1  equations  simply  by  applying  the  Cramer’s 
rule,  it  follows  that 


Po,q 

0 

0 

0 

1 

Po,q 

0 

0 

AV,qix) 

P2,q 

0,  A. 

f3o,q 

0 

^2,qix) 

f^n—l,q 

[  1  ]g  f^n-2,q 

[  2  ]q  (^ri-3,q 

f3o,q 

^n-l,qix 

(^n,q 

[l]q  (^n-l,q 

[2]^ 

An,qix) 

Now,  bringing  the  (n  +  l)-th  column  to  the  first  place  by  n  transpositions  of 
adjacent  columns  and  in  view  of  the  fact  that  the  determinant  of  a  square  matrix 
is  the  same  as  that  of  its  transpose,  second  part  of  assertion  (13)  is  proved.  ◄ 


4.2.  Determinant  forms  for  the  2-iterated  and  hybrid  g-Appell 
polynomials 

We  know  that  the  Bernoulli  polynomials  Bn,q{x)  and  Euler  polynomials 
En,q{x)  are  the  two  particular  members  of  the  g-Appell  family  An^q{x).  First,  we 
establish  determinant  forms  for  the  polynomials  Bn,q{x)  and  En,q{x)  by  choosing 
suitable  values  of  the  coefficients  /3o,q  and  jSi^q  (i  =  1,  2,  •  •  •  ,  n)  in  the  determinant 
form  of  An^q{x). 

Taking  (3o^q  =  1,  (3i^q  =  (^  =  !>  2,  •  •  •  ,  n)  in  the  determinant  form  of  the 

(7-Appell  polynomials  An^q{x)  [5,  p. 12(19)],  the  following  determinant  form  of  the 
(7-Bernoulli  polynomials  Bn,q{x)  is  obtained: 
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Definition  1. 

The  q-Bernoulli  polynomials  Bn,q\ 

(x)  of  degree 

n  are 

to 

o 

II 

=  1, 

En,qix)  = 

= 

1 

X 

1 

x"' 

1 

1 

Wq 

Mg 

1 

Rg 

1 

\n+l\q 

0 

0 

1 

^  •  • 

rnl  1 

0 

=  (-1)" 

LlJg 

1  1  Jq  [n—l\q 

[I 

5 

0 

0 

1 

■  1 

rn-i  1 

1  2  Jg  [n—2]q 

L2jg  [n—l]q 

0 

0 

0 

1 

r  n  ]  1 

In-llg 

defined  by 


n  =  . 


(17) 

The  partieular  case  of  Definition  1,  for  n  =  A,  has  been  eonsidered  in  [15,  p.250]. 

Next,  taking  /3o,q  =  1,  fii,q  =  ^  (i  =  1,  2,  •  •  •  ,  n)  in  the  determinant  definition 
of  the  (7-Appell  polynomials  An^q{x)  [5,  p. 12(19)],  the  following  determinant  form 
of  the  g-Euler  polynomials  En,q{x)  is  obtained: 


Definition  2.  The  q-Euler  polynomials  En,q{x)  of  degree  n  are  defined  by 


Eo,q{x)  =  1, 


En,q{x)  =  (-1)^ 


1  X 


1 


1 
2 

0  1 
0  0 


1 

2 

1  12- 


-2- 

-ij, 


2  LiJ  q 
1 


«n— 1 


1  rn-ln 

2  1  1  Jg 

1  rn-ln 

2  1  2 


1  M 

2  LlJg 

1  w 

2  L2jg 


,  n  =  l,2. 


0  0  0 


1 


(18) 


Remark  4.  Taking  suitable  values  of  the  coefifieients  /3o,g  and  fii^q  (i  =  1,  2,  •  •  •  ,  n) 
in  the  determinant  form  of  the  2IqAP  family,  the  determinant  forms  for  the  2- 
iterated  q-members  and  hybrid  q-speeial  polynomials  ean  be  obtained. 
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First,  the  determinant  forms  for  the  members  of  the  2IgAP  given  in  Table  2  are 
obtained.  Taking  /3o,q  =  1,  (^  =  1;  2,  •  •  •  ,  n)  and  =  Bn,q{x)  in 

f2l 

equation  (13),  the  following  determinant  form  of  the  2IgBP  Bnjg{x)  is  obtained: 


\2] 

Definition  3.  The  2-iterated  q-Bernoulli  polynomials  Bn;q{x)  of  degree  n  are 
defined  by 


B'Sqix)  = 

1  Bi^q{x)  B2,q{x)  ••• 

1  J_  J_  ... 

^  [2],  [3], 

0  1  \^]  w-  ••• 

_  LlJij  12], 

00  1  ••• 

00  0  ••• 

where  Bn,q{x)  (n  =  0, 1,  2, . . .)  are  the 


Bn^ 

-i,<?(^) 

Bn,q{x) 

1 

1 

[n]q 

[n+l]q 

1 

LlJq  [n]q 

q  [n-l]^ 

r2^] 

1 

q  [n—2\q 

L2jq  [n-l\q 

1 

\  ^  ]  -A- 

Ln-lJq  [2], 

(19) 


q-Bernoulli  polynomials  of  degree  n. 


Taking  /3o,g  =  1,  fii^q  =  ^  (f  =  1,  2,  •  •  •  ,  n)  and  Al/g{x)  =  En,q{x)  in  equation 

\2] 

(13),  the  following  determinant  form  of  the  2IgEP  En,q{x)  is  obtained: 


\2] 

Definition  4.  The  2-iterated  q- Euler  polynomials  En,q{x)  of  degree  n  are  defined 
by 
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E&x)  = 

=  1, 

E^nlix)  = 

1 

El,g{x) 

E2,g{x)  ■ 

En—l,g{E) 

En,g{x) 

1 

1 

1 

1 

1 

1 

0 

2 

2 

2 

2 

=  (-ir 

0 

1 

im 

2  lllg 

2  1  1  Ig 

2  lllg 

0 

,  n  =  1,2 

0 

0 

1 

2  1  2  Ig 

1  m 

2  l2lg 

0 

0 

0 

1 

if  ^  1 

2  In-llg 

where  En^g 

(x) 

(n  =  0, 1 

,  2, . . .)  are 

the  q-Euler  polynomials  of  degree  n. 

Next,  the  determinant  forms  for  the  hybrid  g-special  polynomials  given  in 
Table  3  are  obtained. 

Replacing  the  powers  of  x,  i.e.  x^,  x^,  x^, . . .  ,x"'  by  the  corresponding  poly¬ 
nomials  Bo,g(x),  Bi^g(x),  . . . ,  Bn,g{x)  in  the  r.h.s.  and  replacing  x  by  the  polyno¬ 
mial  Bi^g{x)  and  using  the  relation  EBn,q{x)  =  BEn,q{x)  =  En,g{Bi^q{x)}  in  the 
l.h.s.  of  the  equation  (18),  the  following  determinant  form  of  the  gBEP  BEn,g{x) 
(or  gEBP  EBn,g{x))  is  obtained; 


Definitions.  The  q- Bernoulli- Euler  polynomials  BEn,q{x)  (or  q-Euler-Bernoulli 
polynomials  EBn,q{x))  of  degree  n  are  defined  by 
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EBo^q{x) 

BEo^q{x)  - 

=  1, 

EBn,q(^x') 

BBn^qix) 

= 

1 

El,q{x) 

B2,qix)  ■ 

••  Br 

i-l,qix) 

En,q{^) 

1 

1 

1 

2 

1 

2 

1 

2 

1 

2 

0 

0 

1 

Ml],  • 

1 

2 

rM], 

1  m 

2  LlJ, 

0 

=  (-1)" 

,  n 

0 

0 

1 

1 

2 

rM], 

2  L2J, 

0 

0 

0 

1 

if  ^  1 

2  Ln-lJq 

(21) 

where  Bn,q{x)  (n  =  0, 1,  2, . . .)  are  the  q-Bernoulli  polynomials  of  degree  n. 

Next,  replacing  the  powers  of  x,  i.e.  x^, . . .  jX*^  by  the  corresponding 

polynomials  HQ^g{x),  Hi^g{x),  . . . ,  Hn^g{x)  in  the  r.h.s.  and  replacing  x  by  the 
polynomial  Hi^g{x)  and  using  the  relation  BHn,g{x)  =  HBn,g{x)  =  Bn,g{Hi^g{x)} 
in  the  l.h.s.  of  the  equation  (17),  the  following  determinant  form  of  the  gHBP 
HBn,q{x)  (or  (7BHP  BHn,g{x))  is  obtained: 

Definition  6.  The  q-Hermite- Bernoulli  polynomials  HBn,q{x)  (or  q-Bernoulli- 
Hermite  polynomials  BHn,g{x))  of  degree  n  are  defined  by 


BHo^q{x)  =  HBo^q{x)  =  1, 

BByi^qix)  —  H  Bn^gi^x')  — 


1 

Hfix) 

H2{x)  ■ 

Hr, 

.-1(2;) 

H 

n(^) 

1 

1 

1 

1 

1 

X 

[2], 

[3], 

[n]q 

Jn 

'+lb 

l)n 

0 

1 

■  M 

1 

q  [n-1]. 

ra 

1 

,  [n]q 

0 

0 

1 

1 

q  [n-2\q 

[2], 

1 

N-1I9 

0 

0 

0 

1 

L-i 

-]<?  w. 

(22) 


where  Hn,g{x)  (n  =  0, 1,  2, . . .)  are  the  q-Hermite  polynomials  of  degree  n. 


A  Numerical  Computation  of  Zeros 


43 


Further,  replacing  the  powers  of  x,  i.e.  x^, . . . ,  x""  by  the  correspond¬ 

ing  polynomials  HQ^g{x),  Hi^q{x),  . . . ,  Hn,q{x)  in  the  r.h.s.  and  replacing  x  by  the 
polynomial  Hi^q{x)  and  using  the  relation  EHn,q{x)  =  HEn,q{x)  =  En^q{Hi^q{x)} 
in  the  l.h.s.  of  the  equation  (18),  the  following  determinant  form  of  the  gHEP 
HEn^q{x)  (or  (7EHP  EHn,q{x))  is  obtained: 


Definition  7.  The  q-Hermite- Euler  polynomials  HEn,q{x)  (or  q- Euler- Hermite 
polynomials  EHn,q{x))  of  degree  n  are  defined  by 


EHo^q{x)  =  HEo^q{x)  =  1, 

EHn,q{x)  =  HEn,q{x)  = 

1  Hi^q{x)  H2,q{x) 


=  (-ir 


0  0 


2  LlJq 

1 


Hn—l,q{ 

'x) 

Hn,q{x 

1 

1 

2 

2 

1  rn-r 

1  rn] 

2  L  1  . 

<? 

2  LlJq 

1  rn-r 

1  rn] 

2  1  2  . 

9 

2  121  q 

1 

if  ^  1 
2  Ln-lJ, 

n  =  1,2,--- 


where  Hn,q{x)  (n  =  0, 1,  2, . . .)  are  the  q-Hermite  polynomials  of  degree  n. 


(23) 


Remark  5.  Taking  x  =  0  in  determinant  definitions  (17)-(18)  of  Bn,q{x)  and 
En,q{x)  and  on  expanding  the  determinants  with  respect  to  first  row  and  using 
suitable  notations  from  Table  1  (I-II),  the  determinant  definitions  of  the  related 
numbers  Bn^q  and  En,q  can  be  obtained. 


\2\ 

Remark  6.  Taking  x  =  0  in  determinant  definitions  (19) -(23)  of  Bii'qix), 

Ek,q{x),  BEn,q{x)  (or  EBn,q{x)),  BHn,q{x)  (or  HBn,q{x))  and  EHn,q{x)  (or 

HEn,q{x) )  and  then  using  suitable  notations  from  Tables  1,  4  and  5  (I-III),  the  de- 

[2l  [2l 

terminant  definitions  of  the  related  numbers  Bn^q,  En'q,  BEn,q  (or  EBn,q),  BHn,q 
(or  H Bn, q)  and  EHn,q  (or  HEn,q)  Can  he  obtained. 
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4.3.  Identities  for  the  2-iterated  g-Appell  polynomials 

f2l 

In  order  to  establish  the  identities  for  the  21(7 AP  An;q{x),  the  following  iden¬ 
tities  for  the  g-Appell  polynomials  are  considered  [5]. 


^n,q{x) 


X 


n 


“  ~  (^n-k,q  Ak^q{x)j, 

n 

~  Xy  \_}^gP'n-—k,q^k,q{x),  n=l,2,-'' 

fc=0 


n  =  1,2,--- 


Replacing  the  powers  of  x,  i.e.  x^  and  x”  by  the  corresponding  polynomi¬ 
als  Ai^q{x)  and  An^q{x)  in  above  equations  and  then  using  equation  (9)  in  the 

f2l 

resultant  equations,  the  following  identities  for  the  2I(7AP  Anjq{x)  are  obtained: 


The  above  examples  illustrate  that  the  operational  correspondence  established 
in  this  article  can  be  applied  to  derive  the  results  for  the  newly  introduced  q- 
special  polynomials  given  in  Tables  2  and  3  from  the  results  of  the  corresponding 
member  belonging  to  the  (7-Appell  family. 
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